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EXPERIMENTAL INVESTIGATION OF HOT-GAS SIDE HEAT-TRANSFER 
RATES FOR A HYDROGEN -OXYGEN ROCKET , 
by Ralph L. Schacht, Richard J. Quentmeyer, and William L. Jones 

Lewis Research Center 

SUMMARY 

An experimental investigation was conducted at the Lewis Research Center to de- 
termine the hot-gas side heat-transfer rates in a rocket nozzle, Transient temperature 
measurements were made at five axial locations in a copper heat-sink nozzle having ex- 
pansion and contraction area ratios of 4. 64 with gaseous hydrogen and liquid oxygen as 
propellants. The nozzle was operated over a range of chamber pressures from 150 to 
1000 pounds per square inch absolute at fuel percentages of 11 to 17, with most measure- 
ments made at approximately 15 percent hydrogen by weight. 

The data are presented in a correlation equation of the form St*Pr*°' ^ = CRe^'®' 
where St*, Pr*, and Re^ are reference Stanton, Prandtl, and Reynolds numbers, re- 
spectively, and C is a constant. Over the Reynolds number range investigated, the re- 
sults showed that no one constant C correlates the data for all axial locations; how- 
ever, the data did correlate well when each station was handled separately. The values 
of C varies from a high of 0. 026 in the chamber to a low of 0. 015 at the throat. This 
throat value is 42 percent lower than the widely used value of C = 0. 026. 

INTRODUCTION 

In the design of convectively cooled nozzles for high-performance chemical and nu- 
clear rocket systems the high heat fluxes and associated high wall temperatures require 
accurate knowledge of the heat flow across the gas boundary. Present design procedures 
involve gas-side heat-transfer correlations which have not been well substantiated by 
rocket heat-transfer data. Rocket-nozzle heat-transfer investigations are reported in 
references 1 to 5. Although the data of references 1 to 5 appear to correlate moderately 
well, the chamber pressure range, propellant combination, nozzle geometry, and in- 
strumentation technique limit the usefulness of the data. 





It was therefore deemed desirable to conduct a gas-side heat-transfer experiment 
over a wide range of chamber pressures with emphasis on the accurate determination of 
gas-side heat-transfer rates in the regime of high heat flux. Although this program is 
directed toward the solution of problems associated with both the cooling side and the 
hot-gas side, the investigation reported herein had the objective of providing useful hot- 
gas-side design information. 

The rocket nozzle used for this investigation had a 5- inch-diameter throat with a 
4. 64 expansion and contraction area ratio. The propellants used were gaseous hydrogen 
and liquid oxygen. The nozzle wall was constructed of copper, which provided heat-sink 
cooling for the transiently conducted tests. Local heat-transfer rates were determined 
at five axial stations through the nozzle by transient temperature measurements using 
instrumented copper rods. The range of chamber pressures covered was approximately 
150 to 1000 pounds per square inch absolute. The range of mixture ratios investigated 
corresponded to 11 to 17 percent hydrogen by weight, with the bulk of the runs being 
made at 15 percent hydrogen. 


APPARATUS AND EXPERIMENTAL PROCEDURE 
Rocket Design 

A copper heat-sink, solid-wall, hydrogen- oxygen rocket was used in the experiments 
to obtain transient temperature data (fig. 1). The nozzle had a 4. 64 contraction and ex- 


r Gaseous hydrogen 



Figure 1. - Copper heat-sink rocket nozzle with coaxial injector. Dimensions given in inches. 
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(b) Injector with copper faceplate. 
Figure 2. - Coaxial 234 element injector. 
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pansion area ratio. The length of the combustion chamber from the injector to the throat 
was 14. 5 inches. The characteristic length L* was 54 inches. The throat diameter 
was 5 inches. The rocket developed 26 000 pounds of thrust at a chamber pressure of 
900 pounds per square inch absolute. 

A coaxial injector (figs. 1 and 2) was used with liquid oxygen and gaseous hydrogen 
as the propellants. The injector had 234 injector elements uniformly spaced in a dish- 
shaped faceplate (2. 6 elements per square inch). Both chamber and injector were de- 
signed with the idea of obtaining stable engine operation and as uniform a distribution of 
mass flux and temperature as possible to eliminate striation effects. 


Instrumentation 

A small copper plug or rod was devised that, when inserted in the wall of the un- 


x 



pure copper rod 


Instrumentation station 1 



CD-8029 


Figure 3. - Instrumentation locations and sealing details. 
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cooled rocket engine, approximated the model of a one-dimensional semi-infinite slab. 

As shown in figure 3, holes were machined normal to the inside wall of the copper nozzle 
at each of the seven measuring stations (five axial locations and two additional circum- 
ferential locations at the throat). Rods of pure, oxygen-free copper with diameters of 
0. 25 inch were used. The rods were threaded at the end for insertion in the machined 
holes. The threaded portion was machined so that its cross-sectional area closely ap- 
proximated that of the 0. 25-inch-diameter rod. This threaded portion was oxidized be- 
fore insertion to minimize conduction between the rod and the chamber. The length of 
the rods was 2. 5 inches, and the thickness of the chamber wall 2 inches. This length 
was picked to minimize the temperature differential between the heat-transfer surface of 
the rod and the surrounding chamber wall because of the larger volume of a radial ele- 
ment compared to the volume of the rod. 

Figure 3 also shows the thermocouple stations on the copper rod. Six Chromel- 
Alumel thermocouples were used. The Chromel and Alumel wires, 0. 008 inch in diam- 
eter, were peened to the surface. Figure 3 shows a detail of the sealing method used at 
the end of the rod. The cavity surrounding the rod was pressurized with room- 
temperature nitrogen to a pressure 5 to 10 percent below that of the local static pressure 
at each measuring station to prevent hot gas leakage past the threaded portion of the rod. 
The pressures were monitored during each run to ensure that no film cooling was ob- 
tained from the pressurizing gas. Cavity pressure remained constant with time after the 
start of a run. 


Data Recording 

Propellant flow, chamber pressure, and rod temperatures were recorded in a digi- 
tized form on a magnetic tape for direct entry into a digital computer. The digital re- 
cording system had basic sampling rates of 4000 words per second. Figure 4 (p. 6) 
shows the basic block of data as taken by the digitizer. Most data parameters were ar- 
ranged in the block so that, with a minimum of five blocks of data, 60- cycle noise could 
be defined. Then a smoothing technique was used to eliminate the 60- cycle noise and 
greatly diminish any random noise. The calculating procedure fit a smooth curve through 
15 readings of the data parameter, and then one report, to be used in the terminal calcu- 
lations, was made for all parameters at a common time. This reduced the amount of 
terminal calculations. Chamber pressure, which was sampled nine times in each block, 
was used as a triggering device for starting and stopping the calculations. 
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First subscript denotes station; second subscript denotes X/y location. 
Figure 4 - Digitizer data block. 


Test Procedure 

The rocket engine (fig. 5) was installed on a test stand located at the Lewis Plum 
Brook Facility. Propellant valves for controlling gaseous hydrogen and liquid oxygen 
were prepositioned before the run and opened to these fixed positions during the run to 
provide the desired values for chamber pressure and percent of hydrogen. Chamber 
pressures were varied from 150 to 1000 pounds per square inch absolute at 11 to 17 per- 
cent hydrogen. The timing and sequencing of the valves were adjusted so that full cham- 
ber pressure was obtained in 0. 02 to 0. 06 second for a typical run as shown in figure 6. 
This step function in chamber pressure or driving temperature allowed a simpler mathe- 
matical model to be used for obtaining heat-transfer coefficients. 


CALCULATION PROCEDURE 
Constant h Method 

The differential equation for one-dimensional heat conduction in a homogeneous 
medium is 
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Chamber pressure, 
Ib/sq in. abs 
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( 1 ) 


k^ = P cH 
ax 2 dt 


with the boundary conditions 


' k ^= h < T aw- T > 


at X = 0 


( 2 ) 


k— = 0 at X = (3) 

ax 


The initial condition is T = T Q at t = 0. (Symbols are defined in the appendix. ) The 
assumption that the heat-transfer rates are proportional to the difference between the 
temperature at the wall and the temperature of the surrounding medium (eq. (2)) seems 
to be justified because the heat transfer occurs primarily by convection, rather than 
conduction or radiation. Radiant heat transfer is estimated to be an order of magnitude 
less than that due to convection for this experiment. It is therefore neglected in the 
theoretical treatment. The measurements, however, include both effects. The solution 
to equation (1) is given in closed form in reference 6. 

Reference 7 has put this solution in the form 


where 


(?(N x ,t*) = erfc 





erfc z = 1 - — 

Jr 



( 4 ) 


while the dimensionless variables are 
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T - T, 


<P = 


T aw ' T 0 


N, 


hX 

k 


t* = — 

kpc 

Equation (4) has the factor 2 added, which was erroneously omitted in reference 7. 

Equation (4) holds until such time as the temperature starts to rise at the position 
^ = ^he thickness of the slab). For shorter times the temperature distribution given 
by equation (4) is approximately represented by that in a semi-infinite solid, the accuracy 
increasing as X = %<£ approaches the hot boundary (£ -► 0). 

When the temperature starts to rise at X = if, equation (4) fails because of boundary 
condition (eq. (3)). Since longer running times are needed, this failure of the solution 
can be remedied by allowing, in addition to the heat source of strength h(T - T) at 
X = 0, an additional heat source of this strength at X = 2 if. Reference 7 shows that, 
with two identical sources located symmetrically about X = if, the conduction of heat in 
the region 0 < X < <£ is governed by the same differential equation and initial boundary 
conditions as in the slab, but the solution given applies during the time when the temper- 
ature at X = 0 is not influenced by the image source at X = 2if. The temperature of 
the slab is therefore accurately represented within this extended time period by the 
superimposed effect of the two sources: 

^( N X» T ) = ^( N X=|if’ T *) + ^[ N x=(2 -$<£' T *] ( 5 ) 

Equation (5) was programed on a digital computer so that an iteration process found 
the h that satisfied the measured <p for a given location £ and time r*. The last 
term, <p[ N x=(2-Z)&’ T *]> was also reported so that cutoff of data could be made when the 
temperature rise at X = 0 was not yet influenced by the image source at X = 2 if. Since 
the engine runs were all necessarily of short duration, equation (5) was sufficient for 
reduction of the data. 

The properties of materials (k, p, and c) may vary appreciably with temperature, 
and a problem arises concerning which temperature in the material should be used in the 
evaluation. Reference 8 shows that the material properties for simple metals for this 
type of computation should be evaluated at a temperature that is about one -fourth of the 
maximum if the initial temperature is zero, or at (T^ =Q - T Q )/4 + T Q if the initial tem- 
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Figure 7. - Proportionality constant and nondimensional 
temperature as functions of time. Distance normal to 
heated surface, 0. 125 inch; k/pc, 0.001151 square foot 
per second. 


perature is nonzero. The tempera- 
ture of the rod at 4 = 0 was approxi- 
mated with finite differences by the 
equation 

T 4=0 = ^ T 4=0.05" 2T 4=0. l + 3 T 4=0. 2 

where T^_q q 5 , T ^ = o. l> and T 4=0. 2 
are measured values. 

Figure 7 shows the results of a 
study that relates the error in a non- 
dimensional temperature <p to the 
percent error in the nondimensional 
heat transfer hX/k. The nondimen- 
sional temperature cp is plotted 
against time for various hX/k val- 
ues. The proportionality constant 
G from the equation 

« = G d^ 
hX/k 

is plotted against time for the same 
hX/k values. Since the limit of 
running time for a heat- sink engine 
is a constant value of cp (T w close 
to melting at the surface), it can be 


seen from figure 7 that G remains essentially constant with hX/k. Thus, even though 
the length of running time goes down as chamber pressure is increased (hX/k is in- 
creased), the accuracy of obtaining h remains the same for a given measurement in- 
accuracy d<p and a constant cp. If the accuracy in obtaining cp is estimated to be 
±5 percent in the range of cp = 0. 15, the percent error in hX/k is ±8 (0. 15) (0. 05) (100) 


or ±6. 0 percent, as shown by arrow in figure 7. 


Integration Method 

An integration method was also used to check the results from the constant h method 
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described. If the slab is considered to be insulated on all surfaces except the heat- 
transfer surface, the time rate of change of the integrated heat content must be equal to 
the rate of heat flow through the surfaces: 


q = - — f y A X Q X dX = - — f* (pAcT) x dX 
A St J Q A St ^ A 


( 6 ) 


In equation (6) the heat content Q x is above some arbitrary reference level. If the slab 
cross-sectional area is constant over its length, the area term does not appear in equa- 
tion (6). If a short time increment after the start of steady-state burning is considered, 
the time rate of change of integrated heat content at the arithmetic average of these two 
times is equal to the difference in heat content divided by the time interval. The heat- 
transfer coefficient can be obtained from the following equation: 


1 0. 5(t J+tg) (rp 


- T ) 
aw w ; 


0. 5(tj+t2) 




In this equation the pc product is evaluated at the measured temperature at X. 


Combustion Temperature 

Before h can be computed by the methods described previously, the driving tem- 
perature (T ) must be determined throughout the nozzle; to determine T the com- 
bustion temperature must be known. Combustion temperature is a function of combustion 
chamber pressure, percent fuel, and combustion efficiency. 

Numerous programs are presently in existence for the calculation of equilibrium 
compositions and other thermodynamic properties of complex chemical systems. The 
programs of reference 9 for thermodynamic properties and of reference 10 for transport 
properties were modified and simplified for a gaseous -hydrogen - liquid-oxygen system 
by Frederic N. Goldberg of the NASA Lewis Research Center. The program was further 
modified to account for combustion efficiency. The procedure used to account for com- 
bustion efficiency and arrive at the proper combustion temperature is described in the 
paragraphs that follow. 

Under the conditions of one -dimensional isentropic equilibrium flow, the specifica- 
tion of chamber enthalpy H c (for a given mixture ratio), combustion chamber pressure 
P c , and throat area per unit weight flow (A/w)^ may lead to an incompatability. If the 
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calculated value of (A/w) A corresponding to the minimum value of the function 
A/w = f(P, S J (where S„ is the entropy of the chamber gases determined by H and 
P c ) differs from (A/w)^, then one of the input specifications must be relaxed. 

Since the measurements considered to be most reliable were those of combustion 
pressure and geometric throat area, these were considered to be invariant. When the 
calculated value (A/w) m ^ n was greater than (A/ w)^, an iteration was used to determine 
a chamber enthalpy, or equivalently a combustion temperature T c , compatible with the 
measured weight flow and the requirement that choking occurs at the geometric throat 
according to the following algorithm: 





S c +1 = S ( T c +1 > 


where the K = 0 conditions correspond to those of specified input enthalpy. Calculations 

were terminated when there was sufficiently close agreement between the calculated and 

measured throat area per unit weight flow. The square root of the ratio of the final com- 

K=0 

bustion temperature to T £ was used as a measure of combustion efficiency. 

If the initial calculated (A/ w) m ^ n was less than that measured, the enthalpy was left 
invariant. (Note that in this case the preceding method would indicate an efficiency 
greater than 100 percent. ) A new weight flow was then calculated with the combustion 
efficiency assumed to be 100 percent: 


w 

new 


(") 

w'th 


(-) 

W min 


w measured 


The average combustion efficiency obtained from all the test data was 98 percent of 
theoretical equilibrium. 
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Transport Properties 


Eckert and Drake, in a summary discussing turbulent boundary layers (ref. 11), 
recommend that the transport properties appearing in the correlation equations be in- 
troduced at a reference enthalpy given by the equation 

H* = H s + 0. 5(H w - H g ) + 0. 22(Pr*) 1 / 3 (H tQt - H g ) 

instead of by temperature when the range of temperature is such that the specific heat 
varies considerably within the boundary layer. Therefore, the nondimens ional heat- 
transfer parameters were computed by introducing the transport properties as a function 
of H* and P g for the hydrogen- oxygen system described. 


s. 

E 






Figure 8. - Typical profiles of wall temperature, heat flow 
rate, and heat-transfer coefficient plotted against area 
ratio. Chamber pressure, 967 pounds per square inch 
absolute; fuel mixture, 16.5 percent hydrogen; time, 
0.462 second. 


RESULTS AND DISCUSSION 

Figure 8 is a plot of the wall 
temperature, heat flux, and heat- 
transfer coefficient for a typical run, 
for which the chamber pressure was 
967 pounds per square inch absolute, 
the hydrogen concentration was 
16. 5 percent, and the time 0. 462 
second. Figure 8 shows that the heat 
flux was indeed in the range of anti- 
cipated nuclear-rocket nozzles, that 
is, q = 20 Btu per square inch per 
second, and T w = 2000° R. 

The results presented in this 
report are all from the constant h 
solution. As running progressed, 
thermocouples at X/if = 0. 6 and 
X/if = 1. 0 were lost and the inte- 
gration method could not be applied. 
For some of the earlier tests, how- 
ever, results from the integration 
method agreed with the constant h 
solution and gave additional support 
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Figure 9. - Trend comparison of theoretical and measured temperature transients in semi- 
infinite solid. (x * 0.125 in.) 


to this method. For a further test of the apparatus and technique a second nozzle fabri- 
cated and instrumented in the same manner was run at two of the chamber conditions. 
Data reproducibility was shown to be good; that is, it fell within the spread of the data 
and thus indicated nothing unique about one particular apparatus. 

Before the results of any experiment can be used, the validity of the results must be 
established. Figure 9 can be used to consider the results. The solid lines are the theo- 
retical lines from equation (4) rewritten as follows to include the pertinent variables 



for hX/k = 0. 05, 0. 075, 0. 1, and 0. 125 and X = 0. 125 inch. Experimental data from 
one thermocouple for the throat station are shown for three chamber pressures. The 
experimental points follow the contours of the theoretical lines fairly well; thus the as- 
sumption of a constant h solution and a one -dimensional model seems justified. 

The data for the range of chamber pressures from 150 to 1000 pounds per square 
inch absolute are presented in figure 10 for each measuring station of the rocket nozzle 
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(a) Correlation equation, St*Pr*°- 7 = CRe2‘®- 2 . 



Station 3 represents the average values 

for the three measuring stations at the 

throat (3, 3A, and 3B). The nondimen- 

♦ *0 7 

sional parameter St Pr * is plotted 

3|e 

against Re^. All transport properties 
were introduced at the reference enthalpy 
and static pressure for each station. 
Chemical equilibrium conditions were 
used. For these plots duplicate data 
points represent calculations made for 
each of the last two common time reports 
for each run. The last two time reports 
give the greatest T - T Q input into the 
parameter <p = (T - T Q )/(T aw - T Q ) and 
therefore have the least error, as shown 
by figure 7 (p. 10). Figure 9 shows that 
h is fairly constant with time. 

In each part of figure 10 two least- 
square fits for the data are plotted, one 
where the constant and the exponent of 
the Reynolds number were both allowed 

to vary in the standard equation of the form St*Pr*°‘ 7 = CRe^ b , the other where only 
the constant is allowed to vary and the exponent b is held at a fixed value of -0. 2. The 
constants and the exponents are shown for each curve. The plots also have the line rep- 
resented by the equation St*Pr*°- 7 = 0. 026 Re*j"°- 2 , which is similar to the simplified 
Bartz equation (ref. 12). 

In order to give a better statistical interpretation of the data spread around the much 
.*-.* 0. 7 ™_*- 0.2 



(b) Correlation equation, $t*Pr*°' 7 - CRe^. 

Figure 11. - Exponent and correlation constant as functions of 
area ratio for two correlation equations. 


used St Pr 


= CRe: 


correlation, the standard deviation and the percent deviation 


for each station are presented in the preceding table. The column headed 1. 96 a shows 
the deviation around C defining the 95 percent confidence interval. 
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Over the Reynolds number range covered by these experiments either correlation 
(St*Pr*®' ^ = CRe*j”^' 2 or St*Pr*^' ^ = CRe*j b ) seems to be adequate for obtaining heat- 
transfer coefficients at a given station for the hydrogen-oxygen system used. 

Figure 11(a) is a plot of the constant against area ratio for an exponent of -0. 2. 

These data indicate that no one value of C will correlate the data for all axial locations. 
In fact, the constant C = 0. 026, which has been widely used to predict gas-side heat- 
transfer coefficients for nozzles, overpredicts the coefficients for all axial locations 
except the chamber. The measured heat-transfer coefficients at the throat were a PP^°^‘ 
imately 42 percent less than those predicted by the equation St*Pr* * = 0. 026 Re d ‘ . 

Figure 11(b) shows plots of the exponent and the constant against area ratio where the 

curve fit of the data allowed both to vary. 

From hydrogen-oxygen data presented in reference 5 for a chamber pressure of 
150 pounds per square inch absolute and a mixture ratio of 5.0 (16. 7 percent hydrogen), 
it was possible to compute values of the throat heat-transfer coefficient h. This com- 
putation gave a band of h from 0. 0014 to 0. 0010 Btu per square inch per second per °R, 
which compares with a value of 0. 00136 measured at the throat of the nozzle tested 
herein under the same operating conditions. Thus, the level of the heat-transfer coef- 
ficients at the throat shows reasonable agreement for the hydrogen- oxygen system at low 
chamber pressures. 

Data taken from reference 3 for three nozzle configurations were recalculated using 
a reference enthalpy input of all transport properties and fitted to the equation 
St*Pr*°‘ ^ = CRe d ”^' 2 . The resulting C values at the throat were 0. 019, 0. 017, and 
0. 023, respectively. Data taken from reference 4 were recalculated in a similar man- 
ner, and a value of C =* 0. 018 at the throat resulted. The curves of C against area 
ratio A/A* also showed trends similar to those of figure 11 (i. e. , values of C high in 
the chamber, low in the throat, and increasing again at the exit). A mention of this work 
is given only to say that, while the variation of C with A/A* that was found is con- 
sidered correct for the nozzle geometry, propellant system, and injectors used, the 
values should not be used universally for calculating the heat-transfer rate to other 
nozzles. 


CONCLUDING REMARKS 

If a designer desires to use the correlations presented herein for a system similar 
to the rocket of this experiment, he must work with the thermodynamic and transport 
properties for chemical equilibrium as computed herein (refs. 9 and 10). All properties 
must be inserted at the reference enthalpy and static pressure of the axial station where 
the heat-transfer coefficient is desired. The effect of using different properties can be 
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great. Using the new transport properties (ref. 10) reduced the correlation constant C 
by approximately 30 percent over those obtained by using transport properties from 
reference 13. Transport properties given in reference 10 are considered to be more 
accurate than those previously used. 

The values of the correlation constant C reported herein are for a given nozzle 
geometry, propellant combination, and injector configuration and are not to be con- 
sidered universal for rocket heat-transfer calculations. 

SUMMARY OF RESULTS 

An experimental investigation was conducted at the Lewis Research Center to de- 
termine the hot-gas-side heat-transfer rates in a rocket nozzle. Transient temperature 
measurements were made at five axial locations in a copper heat-sink nozzle with ex- 
pansion and contraction ratios of 4. 64 with gaseous hydrogen and liquid oxygen as pro- 
pellants. The nozzle was operated over a range of chamber pressures from 150 to 
1000 pounds per square inch absolute at fuel percentages ranging from 11 to 17, with 
most measurements made at approximately 15 percent hydrogen. The following results 
were obtained: 

1. The data at a given location in the nozzle are adequately correlated by the equa- 

3|c 5|c0 n *-0 2 

tion St Pr ' = CRe^ ' , where St , Pr , and Re^ are reference Stanton, Prandtl, 

and Reynolds numbers, respectively, and C is a constant. 

2. The constant C varies with axial location and is less than 0. 026 for all locations 
except the chamber. 

3. The constant C at the throat was 42 percent less than the widely used value of 
C = 0. 026. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 16, 1965. 
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APPENDIX - SYMBOLS 


A cross-sectional area, sq in. 

b exponent in equation 

St*Pr*°- 7 = CReJ b 

C constant in equation 

St*Pr* 0 ' 1 ■= CRe^ 

c specific heat of materials, 

Btu/(lb)(°R) 

G proportionality constant; 



St Stanton number 

T temperature, °R 

T a» { < H aw P s> 

t time, sec 

tj time 1, sec 

tg time 2, sec 

V volume, cu in. 

w total weight flow rate, lb/sec 

w 0 oxygen weight flow rate, lb/sec 

°2 


H enthalpy 

H aw adiabatic wall enthalpy, 

H s + (H tot - H s ) 

h heat-transfer coefficient, 

Btu/(sec)(sq in. )(°R) 

k thermal coefficient of conduc- 

tivity of material, 

Btu/(sec)(in. )(°R) 

L* characteristic length of rocket 

engine, in., Vg/A^ 

length of heat flow path in rod, in. 

N x nondimensional parameter, hX/k 

P pressure, lb/sq in. abs 

Pr Prandtl number 

Q quantity of heat, Btu/ cu in. 

q heat flow rate per unit area, 

Btu/(sec)(sq in. ) 

Re^ Reynolds number based on di- 

ameter 

S entropy, Btu/(lb)(°R) 


X distance normal to heat surface, 

in. 

£ nondimensional parameter, X/fif 

p material density, lb/cu in. 

a standard deviation 

r* nondimensional parameter, 

h 2 t/kpc 

(p nondimensional temperature, 

(T - T 0 )/(T aw - T 0 ) 

^ average value for three circum- 

ferential measuring stations 
at throat of rocket (station 3) 

Subscripts: 

A, B circumferential locations at 

throat (station 3, 120° apart) 

c chamber 

min minimum 

orf orifice 

s static 

th throat 
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tot 

total 

12 3 

’ ’ instrumentation stations 

w 

wall 

4, 5 

X 

distance 

Superscript: 

0 

zero burning time or start of 
burning 

* reference condition 
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